Matter Bounce Cosmology with the /(T) Gravity 
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We show that the /(T) gravitational paradigm, in which gravity is described by an arbitrary 
function of the torsion scalar, can provide a mechanism for realizing bouncing cosmologies, thereby 
avoiding the Big Bang singularity. After constructing the simplest version of an /(T) matter bounce, 
we investigate the scalar and tensor modes of cosmological perturbations. Our results show that 
metric perturbations in the scalar sector lead to a background-dependent sound speed, which is a 
distinguishable feature from Einstein gravity. Additionally, we obtain a scale-invariant primordial 
power spectrum, which is consistent with cosmological observations, but suffers from the problem 
of a large tensor-to-scalar ratio. However, this can be avoided by introducing extra fields, such as a 
matter bounce curvaton. 
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I. INTRODUCTION 

Inflation is now considered to be a crucial part of 
the cosmological history of the universe [l[ , however the 
so called "standard model of the universe" still faces 
the problem of the initial singularity. Such a singular- 
ity is unavoidable if inflation is realized using a scalar 
field while the background spacetime is described by the 
standard Einstein action [2[- As a consequence, there 
has been much energy expended towards resolving this 
problem, e.g. through null-energy-condition violating 
quantum fluctuations ( "island cosmology models" d, 0] ) , 
quantum gravity effects, or effective field theory tech- 
niques. 

A potential solution to the cosmological singularity 
problem may be provided by non-singular bouncing cos- 
mologies 0. Such scenarios have been constructed 
through various approaches to modified gravity, such 
as the Pre-Big-Bang and the Ekpyrotic Q models, 
gravitational actions with higher order corrections M, 
braneworld scenarios ^ , non-relativistic gravity [l^, [HI , 
loop quantum cosmology [l2| or in the frame of a closed 
universe [l^l • For a review on models of modifying higher 
derivative gravity, which solve cosmic singularities in an 
efficient way, we refer to Ref. ■ Non-singular bounces 
may be alternatively investigated using effective field the- 
ory techniques, introducin g m atter fields which violate 
the null energy condition Il5|. [l6j or introducing non- 
conventional mixing terms |17| . The extension of all the 
above bouncing scenarios is the (old) paradigm of cyclic 
cosmology (l8| , in which the universe experiences a peri- 
odic sequence of contractions and expansions, which has 



'Electronic address; |ycai21@asu.edu] 
t Electronic address: schcnl02@asu.e dul 
tElcctronic address: jbdent@asu.edu' 
^Electronic address: sourisli.d@gmail.com 
^Electronic address: ,EmmanuerSaridakis@baylor.edu| 



been reawakened recently |i9j . since it brings different 
insights to the origin of the observable universe pol - [23j 
(see [i^l for a review). 

Along separate lines, f(T) gravity has recently received 
attention in the literature [25-[27|, mostly in the context 
of explaining the observed acceleration of the Universe. 
It is based on the old idea of the "teleparallel" equivalent 
of General Relativity (TEGR) IMIli, which, instead of 
using the curvature defined via the Levi-Civita connec- 
tion, uses the Weitzenbock connection that has no cur- 
vature but only torsion. The dynamical objects in such 
a framework arc the four linearly independent vierbeins. 
The advantage of this framework is that the torsion ten- 
sor is formed solely from products of first derivatives of 
the tetrad. As described in the Lagrangian den- 
sity, T, can then be constructed from this torsion tensor 
under the assumptions of invariance under general coor- 
dinate transformations, global Lorentz transformations, 
and the parity operation, along with requiring the La- 
grangian density to be second order in the torsion ten- 
sor. However, instead of using the tensor scalar T the 
authors of [2y, generalized the above formalism to a 
modified f{T) version, thus making the Lagrangian den- 
sity a function of T, similar to the well-known exten- 
sion of f{R) Einstein-Hilbcrt action. In comparison with 
f{R) gravity, whose fourth-order equations may lead to 
pathologies, f{T) gravity has the significant advantage 
of possessing second-order field equations. This feature 
has led to a rapidly increasing interest in the literature, 
and apart from obtaining acceleration (26l . [27j one can 
reconstruct a variety of cosmological evolutions [131 and 
solutions m, add a scalar field [Hj, use observational 
data in order to constrain the model parameters 13311 , ex- 
amine the dynamical behavior of the scenario |34l |. and 
proceed beyond the background evolution, investigating 
the vacuum and matter perturbations (sBI . [36j as well as 
the large-scale structure [s^]- Note however that there 
is a discussion whether forms of /(T) other than linear 
should be expected [ssj . 

One interesting feature of the f{T) theory is that the 
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null energy condition could be effectively violated. Ac- 
companied with this feature, one expects to obtain a sc- 
ries of nontrivial phenomenological solutions. Particu- 
larly, it was observed that the violation of null energy 
condition is related to the nonsingular bouncing solu- 
tion in the early universe [isj . Therefore, the singularity 
avoidance can be obtained in general in /(T) gravity. In 
the present work we are interested in searching for a sce- 
nario of bouncing cosmology in the early universe, which 
is governed by f(T) gravity. As we show, the realiza- 
tion of a big bounce and the avoidance of singularities 
is straightforward. Additionally, we investigate in detail 
the evolution of perturbations in a specific model of the 
matter bounce. Our example illustrates that it is possible 
to generate a scale-invariant primordial power spectrum 
and pass through the nonsingular bouncing point in the 
context of /(T) bounce cosmology. 

This paper is organized as follows. In section |TT] we 
briefly review the basic idea of f{T) gravity. In section 
mil we investigate the realization of bouncing cosmology 
by virtue of /(T) gravity. Specifically, we postulate an 
explicit form for the background scale factor, and follow- 
ing this ansatz we reconstruct the corresponding form of 
f{T) both analytically and numerically. In section IIVI 
we study the evolution of cosmological perturbations of 
scalar and tensor types along with the matter bounce 
scenario. Finally, section |V] is devoted to the summary 
of the results. 



II. f{T) GRAVITY AND COSMOLOGY 

In this section we briefly review /(T) gravity and 
we provide the background cosmological equations in a 
universe governed by such a modified gravitational sec- 
tor. Throughout the work we consider a flat Friedmann- 
Robertson- Walker (FRW) background geometry with 
metric 

ds^ = dt^ - a^{t) S^jdx'dx^ (1) 

where t is the cosmic time, are the comoving spatial co- 
ordinates, and a{t) is the scale factor. In this manuscript 
our notation is as follows: Greek indices n,i>,... run over 
all coordinate space-time 0, 1, 2, 3, lower case Latin in- 
dices (from the middle of the alphabet) i, j, ... run over 
spatial coordinates 1, 2, 3, capital Latin indices A,B,... 
run over the tangent space-time 0, 1, 2, 3, and lower case 
Latin indices (from the beginning of the alphabet) a, b,... 
will run over the tangent space spatial coordinates 1, 2, 
3. 

A. /(T) gravity 

As stated in the Introduction, the dynamical variable 
of the old "teleparallel" gravity, as well as its /(T) ex- 
tension, is the vierbein field eAix'^). This forms an or- 
thonormal basis for the tangent space at each point x'^ 



of the manifold, that is e^i • = tiab, where tjab = 
diag{l, —1, —1, —1). Furthermore, the vector can be 
analyzed with the use of its components in a coordi- 
nate basis, that is = e^9^. 

In such a construction, the metric tensor is obtained 
from the dual vierbein as 

g^,ix)=7^ABef,ix)e^(x). (2) 

Contrary to General Relativity, which uses the torsion- 
less Levi-Civita connection, in the present formalism 
ones uses the curvatureless Weitzenbock connection [s^ , 
whose torsion tensor reads 

= r.,, - = e\ {d^et - d.e^). (3) 

Moreover, the contorsion tensor, which equals the differ- 
ence between the Weitzenbock and Levi-Civita connec- 
tions, is defined as 

K^"'p = -\[t^"p-T''''p-T,^''). (4) 

Finally, it proves useful to define 

S^^^^\ (a-'-; + T^l 5'; T"^i) . (5) 

Using these quantities one can define the so called 
"teleparallel Lagrangian", which is nothing other than 
the torsion scalar, as [H, [23, El 

r^5/^rv- (6) 

In summary, in the present formalism all the information 
concerning the gravitational field is included in the tor- 
sion tensor T^^, and the torsion scalar T arises from it 
in a similar way as the curvature scalar arises from the 
curvature (Riemann) tensor. Finally, the torsion scalar 
gives rise to the dynamical equations for the vierbein, 
which imply the Einstein equations for the metric. 

While in teleparallel gravity the action is constructed 
by the teleparallel Lagrangian T, the idea of f{T) gravity 
is to generalize T to a function T + f{T), which is similar 
in spirit to the generalization of the Ricci scalar R in the 
Einstein-Hilbert action to a function f{R). In particular, 
the action in a universe governed by f(T) gravity reads: 

I = I d'xe [T + f{T) + L„] , (7) 

where e = dct(e^) ~ \/~-g and Lm stands for the matter 
Lagrangian. We mention here that since the Ricci scalar 
R and the torsion scalar T differ only by a total derivative 
[i^ . in the case where f(T) is a constant (which will 
play the role of a cosmological constant) the action ([7]) 
is equivalent to General Relativity with a cosmological 
constant. 

Lastly, we stress that throughout this work we use the 
common choice for the form of the vierbien, namely 

= A\&g{l,a,a,a). (8) 
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It can be easily found that the family of vierbiens related 
to ^ through global Lorentz transformations, lead to 
the same equations of motion. Note however that, as 
was shown in [i^, /(T) gravity does not preserve local 
Lorentz invariance. Thus, one should in principle study 
the cosmological consequences of a more general vierbicn 
ansatz, but for simplicity we remain with choice (jSJ (see 
also dl). 



where /.t and f^TT denote respectively the first and sec- 
ond derivatives of the function f{T) with respect to T, 
and the mixed indices are used as in Sa'^'^ = e'^^Sp'^". 

em 

Note that the tensor T on the right-hand side is the 
usual energy-momentum tensor. 

If we assume the background to be a perfect fluid, 

em 

then the energy momentum tensor takes the form T py = 
V9tiv — [p + p)upUy, where is the fluid four- velocity. 
Under this, one sees that equations © lead to the back- 
ground (Friedmann) equations 

= —P,n - ^ - 2/,Ti?' (10) 
^ 4:TTG{p,n + P,n) 

1 + /r - 12i/2/,TT ■ ^ ' 

In these expressions we have introduced the Hubble pa- 
rameter H = a/a, where a dot denotes a derivative with 
respect to coordinate time t. Moreover, pm and Pm stand 
respectively for the energy density and pressure of the 
matter content of the universe, with equation-of-state 
parameter Wm = Pm/ Pm- Finally, we have employed the 
very useful relation 

^=-6i^^ (12) 

which straightforwardly arises from evaluation of ([5]) for 
the unperturbed metric. 

III. THE BACKGROUND SOLUTION OF /(T) 
MATTER BOUNCE 

In this section we examine how cosmological scenar- 
ios governed by f(T) gravity can produce a cosmological 
bounce. There are two distinct points in such an investi- 
gation. The first is to examine whether the background 
evolution allows for bouncing solutions. If this is indeed 
possible, then the second point is to examine the evolu- 
tion of perturbations through the bounce. The first task 
is the subject of this section, while the second one will be 
investigated in the next section. Finally, we mention that 



B. Background /(T) cosmology 



Let us now present the background cosmological equa- 
tions in a universe governed by f(T) gravity. Variation 
of the action ([7]) with respect to the vierbcin gives the 
equations of motion 



(9) 

I 

in order to be closer to the convention of the literature 
on this field, we use Mpi instead of G when necessary, 
using the relation Mpi = 1/ VSttG. 

In principle, whether a universe is expanding or con- 
tracting depends on the positivity of the Hubble param- 
eter. In the contracting phase that exists prior to the 
bounce, the Hubble parameter H is negative, while in 
the expanding one that exists after it we have H > 0. 
By making use of the continuity equations it follows that 
at the bounce point H = 0. Finally, it is easy to see that 
throughout this transition H > 0. On the other hand, 
for the transition from expansion to contraction, that is 
for the cosmological turnaround, we have H > before 
and H < after, while exactly at the turnaround point 
we have H ^ 0. Throughout this transition H < 0. 

Having in mind the above general requirements for a 
cosmological bounce, and observing the Friedmann equa- 
tions pO)) . pT|) . we deduce that such a behavior can be 
easily obtained in principle in the context of f(T) cos- 
mology. In particular, one starts with a specific, desir- 
able form of the bouncing scale factor a(t), and thus one 
immediately knows H[t). Concerning the matter fluid 
content of the universe, with equation-of-state parame- 
ter Wmi its evolution equation p„i + ZHil -t- 'Wm)Pm — 
straightforwardly gives the solution Pm{i) since a{t) is 
known. Thus, we can insert these relations in (jlOp . and 
determine /(T), that is f{—6H'^), which generates such 
an H{t) solution. Although this procedure can always 
be performed numerically, in the following we present a 
simple example of a bouncing solution that allows for 
analytical results. 

We start with a bouncing scale factor of the form 

/ 3 \ 
a(t) = as h + -CT<M , (13) 

where is the scale factor at the bouncing point, and 
cr is a positive parameter which describes how fast the 
bounce takes place. Such an ansatz presents the bounc- 
ing behavior, corresponding to matter-dominated con- 
traction and expansion, and additionally it exhibits the 
advantage of allowing for semi-analytic solutions. In such 



e-'dpieSAm^ + /,t] - e^AT'.xSp-^' + S a^" d p{T) J ^tt - -eW + /(T)] = 4^Ge^ Tp^ 
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a scenario t varies between —oo and +00, with t = the 
bounce point. Finally, in the following we normalize the 
bounce scale factor as to unity. 
Straightforwardly we find 



Hit) 



at 



(1 + 3(7^2/2) 



Tit) = - 



1 



^CTt^?' 



Therefore, provided —\Jf^ \J inversion of 

this expression gives the t(r) relation as 



t{T) = ± 



4 

3T 



2 

3^ 



4^/IV 



3Tcr2 



1/2 



(14) 



where we have kept the solution pair that gives the cor- 
rect (t = at T = 0) behavior. Notice that when t > 

\ hr- and t < — 1 / we have assumed the usual Einstein 

y OCT W OCT 

gravity or the TEGR to be the prevailing framework, 
thus negating the need to pursue an / (T) action in that 
region. Furthermore, we assume the matter content of 
the universe to be dust, namely w,„ ~ 0. When inserted 
in the evolution-equation, this leads to the usual dust- 
evolution pm ~ PmBO^glo-^ 1 with p„iB its valuc at the 
bouncing point. 

Inserting the above expressions into (|10p we obtain a 
differential equation for fit), which can be easily solved 
analytically as 



fit) 



it 



(2 + 'Saf^)M^i 



6tM^ 
t 2 + 3t^a 



Go'PmB ArcTan 




(15) 



We mention that in the calculation we have set the inte- 
gration constant to be zero so that the solution is con- 
sistent with the Friedmann equation. Thus, the corre- 
sponding /(T) expression that generates a bouncing scale 
factor of the form (|13p arises from expression (jlSp with 
the insertion of the t(r) relation from (|T4)) . Note that 
the solution (jlSP is an even function of t, and thus the ± 
solutions of (fT4|) correspond to the contraction and ex- 
pansion phase respectively. Obviously they give the same 
form of /(T). 

In order to present this behavior more transparently, 
in Fig. [1] we depict /(T) that generates the matter- 
dominated bounce in /(T) gravity, with = 1, cr = 
7 X 10-^M|., and p„B = 1.41 x 10-^Af|,,. Note that, the 
value of a is roughly determined by the amplitude of the 
CMB spectrum, and that of PmS depends on how fast 
the standard Einstein gravity is recovered in the /(T) 
theory. Their physical meanings will be discussed in the 
next section in detail. 

Furthermore, we numerically derived the evolution of 
the /(T) and the Hubble parameter H as functions of the 
cosmic time in Figs. [2] and [3] respectively. Particularly, 
Fig. [2]shows that the evolution of /(T) is symmetric with 
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FIG. 1: The form of f as a function of the torsion scalar T 
in matter bounce cosmology. The parameters a and pmB were 
chosen to be a ^ 7 x W'^M^i and pmB = 1.41 x 1Q~^ M^i 
respectively, and the graph is in units of Mpi . 
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FIG. 2: Evolution of the function f in terms of the cosmic 
time t in matter bounce cosmology. The parameters a and 
pmB were chosen to be a — 7 x 10~^Mpi and pmB ~ 1.41 x 
W~"'Mpi respectively, and the graph is in units of Mpi. 



respect to the bouncing point t = 0. At the bouncing 
point / arrives at a minimal value WirGprnB, which hap- 
pens to cancel the contribution of normal matter fields, 
and thus leads to the nonsingular bounce. From Fig. [3l 
one can read that the background evolution of the uni- 
verse follows the usual Einstein gravity away from the 
bouncing phase, but it is dominated by /(T) in the mid- 
dle period. These feature are completely consistent with 
the designs of the scenario as expected. 
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FIG. 3: Evolution of the Hubble parameter H in terms of the 
cosmic time t in matter bounce cosmology. The parameters 
a and pmB were chosen to be a = 7 x 10~^Mpi and pmB ~ 
1.41 X 10~^Mpj respectively, and the graph is in units of Mpi. 



IV. COSMOLOGICAL PERTURBATIONS IN 
THE /(T) MATTER BOUNCE 

In the previous section we presented a simple realiza- 
tion of the eosmological bounee in /(T) cosmology, at the 
background level. In this section we extend our analysis 
to an investigation of the perturbations. 

We begin with a brief discussion of the eosmological 
evolution of primordial perturbations in the framework 
of a flat FRW universe. A standard process for gener- 
ating a primordial power spectrum suggests that eosmo- 
logical fluctuations should initially emerge inside a Hub- 
ble radius, then exit it in the primordial epoch, and fi- 
nally re-enter at late times. This process can be achieved 
in the matter bounce cosmology (see for example (45|). 



In this scenario there exist quantum fluctuations around 
the initial vacuum state, well in advance of the time 
when the bouncing point is reached. Along with the 
matter-like contraction, the quantum fluctuations would 
exit the Hubble radius, since the Hubble radius decrease 
faster than the wavelengths of the primordial fluctua- 
tions. When passing through the bouncing point, the 
background evolution could affect the scale dependence 
of the perturbations at ultra-violet scales. However, the 
observable primordial perturbations, responsible for the 
large scale structure of our universe, are mainly origi- 
nated in the infrared regime, where the modified grav- 
ity effect becomes very limited [46j . Consequently, in 
a generic case, one can estimate the formation of pri- 
mordial power spectrum with the standard eosmological 
perturbation theory. In the following, we study the per- 
turbation theory in f(T) cosmology in detail and we ver- 
ify this statement in a specific model of matter bounce 
cosmology. 

A. Generic analysis 

To begin with, we shall work in the longitudinal gauge 
which only involves scalar-type metric fluctuations as 

ds^ = (1 + 2$)dt^ - a^it){l - 2*)(ix^ , (16) 

thus, as usual, the scalar metric fluctuations are charac- 
terized by two functions $ and Correspondingly, the 
fluctuation of torsion scalar at leading order is given by 

ST ^ 12H{^ + H^) , (17) 

which will be widely used in the following calculation. 

By expanding the gravitational equations of motion 
to linear order, we obtain the following perturbation 
equations [11] [s^ . 



V2 

(l + /,T)-^^'-3(l + /,T)i^*-3(l + /,T)i^^$ + 36/,TTi^^(* + i^3') =4nGSp , (18) 

(l + /.T-12i?V,TT)(* + i/$) =4TTGSq , (19) 

(1 + /.t)(*-*) =87rG^s , (20) 

(1 + It - 12i/2/.Tr)* + H{1 + /,t - UH^Itt)^ 

+ 3H{1 + f^T - 12i/V,TT - ^2Hf^TT + ^^H^Hf^TTT)^ 

+ [3iJ^(l -I- /,T - I2i?V,TT) + 27j(l + f^T - 30H^f,TT + 72ijV,TTT)]^> 

+ ^M'^ V^(^ - $) =inG5p . (21) 

2a^ 



The functions Sp, Sp, Sq, and 5s are the fluctuations of stress, respectively. We take the matter component to 
energy density, pressure, fluid velocity, and anisotropic be a canonical scalar field (j) with a Lagrangian in form 
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of 



and thus wc acquire 



5p 
Sq 
Ss 
Sp 



<j>id<p - 0$) + V^5(f> , 



= , 

= 0(50 - 0$) - F^5(/> , 



(22) 
(23) 
(24) 
(25) 



respectively. 

Inserting relation ([24]) into ([20]) one can obtain \1/ = $ 
due to a vanishing anisotropic stress, which is also widely 
found in the standard theory of cosmological perturba- 
tions (for example see [13]). Moreover, combining p9)) 
and ([231) implies that the gravitational potential $ can be 
completely determined by the scalar field fluctuation Scj). 
Therefore we conclude that for our choice of the tetrad 
given in ([8]) , there exists only a single degree of freedom 
in the scenario of f{T) gravity minimally coupled to a 
canonical scalar field. Note that there is another evolu- 
tion equation to describe the dynamics of cosmological 
perturbations, namely the perturbed equation of motion 
for 50. However, using (|19p it can be shown that this 
second equation is consistent with ([T5)) and (PT|) . 

In order to understand the evolution of scalar-sector 
metric perturbations, we use the perturbed equation of 
motion for the gravitational potential $ instead of the 
scalar field fluctuation Sep. Combining ([T^ . (PT|) and 
we obtain the complete form of the equation of motion 
of one Fourier mode $/j as: 



3>fc + a$fe + M-^fc 



c 







with 



.TTT) 



6H^ + 2H 



1 + /,T 

2HV^ 



l2H^f,TT 



i + It- 

1 + /,T 



12i72/ 



TT 



1 + /,T - UH^f^TT 



(26) 
,(27) 

(28) 
(29) 



where we have applied the relation 4* = $ holding in the 
absence of anisotropic stress. The functions a, fj? and 
Cg are respectively the frictional term, the effective mass, 
and the sound speed parameter for the gravitational po- 
tential 4>. Moreover, we recall the scalar-field background 
equation 



3H(I) + V^ = 0, 



(30) 



and the second Friedmann equation (jlip . which in our 
case reads: 



(1 + f,T - l2H^f^TT)H = -47rG02 



(31) 



Consequently, to make use of these two background equa- 
tions, (pS)) can be further simplified as 



HH 




2H-^\^L. + = 0. 



H 

(32) 

Surprisingly, we find that the equation of motion for the 
gravitational potential in the present /(T) scenario is the 
same with that in the standard Einstein gravity [13] , ex- 
cept the newly introduced sound speed parameter. This 
important feature could be a key for us to explore poten- 
tial clues of the f{T) theory in cosmological surveys. 



B. Variables of perturbations in bouncing 
cosmology 

After having derived the equation of motion for the 
gravitational potential, we proceed in solving it in the 
detailed bouncing cosmology. As we showed in the back- 
ground analysis of section IIIIl we can achieve the exact 
matter bounce scenario in which the universe evolves as 
a matter-dominated one in the contracting phase. In or- 
der to accommodate this picture, the matter component 
could be a massive scalar field similar to the mechanism 
of the Lee- Wick bounce [i^ , or a scalar field with a fine- 
tuned exponential potential (48| . 

One often uses a gauge-invariant variable C, the curva- 
ture fluctuation in comoving coordinates, to characterize 
the cosmological inhomogeneities. In the case of /(T) 
cosmology we assume that the form of C is the same as 
that defined in the standard cosmological perturbation 
theory, which is given by 



TT 



(33) 



A useful relation for the time derivative of C can be de- 
rived upon making use of equation p2p . namely 



H 



(34) 



In a generic expanding universe approaches zero at 
large length scales, fc — >■ 0, since the dominant mode of <I>fc 
is then nearly constant. However, in the matter bounce 
cosmology the metric perturbation $fc in the contracting 
phase is dominated by a growing mode with $/; ~ fc"^/^, 
and thus C ke eps increasing before arriving at the bounc- 
ing point |49j . 

Note that, one may be concerned of the variable ( 
becomes ill-defined when H changes sign. In fact, this 
specious trouble was extensively studied in many aspects 
of cosmological perturbation theory. At present, our un- 
derstanding on a well-defined cosmological perturbation 
theory is to require the metric perturbation and the cor- 
responding extrinsic curvature behave smoothly through- 
out the background evolution. The pioneer discussion on 
this topic appeared in Ref. [11] , and we refer to Ref. [1^ 
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for recent detailed analysis by tracking the evolution of 
metric perturbation step by step in matter bounce cos- 
mology. In our calculation we still make use of C merely 
since its analytic analysis is very convenient to be per- 
formed away from the bouncing phase. In addition, we 
would like to point out that the knowledge obtained in 
General Relativity can be also applied to /(T) gravity, 
which is a modification of Einstein gravity, in the frame of 
bouncing cosmology. The reason is that the observable 
modes of perturbations in bouncing cosmology are dis- 
tributed in the Infrared regime and thus the effect caused 
by the modification of Einstein gravity has to be very 
limited. 

In order to simplify the calculation, we further intro- 
duce a canonical variable to characterize the cosmological 
perturbations 



where 



with e — - 
of motion 



aV2e 



(35) 



(36) 



^ It can then be found that the equation 



TP 



-)vk=0 



(37) 



where the prime denotes the derivative with respect to 
the comoving time r = J dt/ a, is still available in a con- 
tracting universe. 



The primordial power spectrum of the /(T) 
matter bounce 



In order to perform a specific analysis we recall that in 
the matter-like contracting phase the scalar factor evolves 



t 



2/3 



(38) 



and z cx a. In this period, by solving the background 
equations of motion, one obtains the approximate rela- 
tions 



-1 



PmB 
2A/V 



(39) 



that hold far before the beginning of the bouncing era 
tm = ~\J'^- As we mentioned above, pmB is the energy 
density of the matter field at the bouncing point and a 
describes how fast the bounce takes place. Therefore, 
the sound speed of the curvature perturbation reverts to 
Cj ~ 1 in the matter-like contracting phase. 

We mention here that inserting the f{T) form of ([M)) 
into the action ([7]), we find that the standard Einstein 
gravity will automatically be recovered when pmB — 



2Mpia. Particularly, when exactly equals to 2Mpia, 
we get f{T) ^ 0{T^), which wiU dilute faster than the 
Ricci scalar during late-time evolution. Even when p^B 
is not equal to 2Mpia, it is clear that the system satisfies 
General Relativity with a rescaled gravitational constant. 
Thus, the combination of pmB and a could, in principle, 
be constrained by measurements of the gravitational con- 
stant. In our computation we choose pmB to be slightly 
different from 2A/pj(T, thus our model is able to approach 
the standard Einstein theory far away from the bounce. 
As a consequence, the perturbation equation becomes 



(40) 



in the contracting phase. Initially the /c^-terni dominates 
in (|40p and thus we can neglect the gravitational term. 
From this point of view, the fluctuation corresponds to a 
free scalar propagating in a flat spacetime, and naturally 
the initial condition takes the form of the Bunch-Davies 
vacuum 15511: 



-ikr 



Vk 



2k 



Making use of the vacuum initial condition we solve the 
perturbation equation exactly, obtaining the solution 



'ikr 



Vk 



1 



kr 



From this result we deduce that the quantum fluctua- 
tions could become classical perturbations, after exiting 
the Hubble radius, due to the gravitational term in equa- 
tion PO)) . Moreover, the amplitude of the metric pertur- 
bations will keep increasing until the universe arrives at 
the bouncing phase at the moment t„i- 

From the definition of the power spectrum we see that 
C ^ k^/^\vk\ is scale-invariant in our model, which can 
also be achieved in infiationary cosmology. However, the 
coefficient e takes the value | in the matter-like contrac- 
tion and thus it is unable to amplify the power spectrum 
of metric perturbation as in inflation. A detailed calcu- 
lation provides the expression of the primordial power 
spectrum for the /(T) matter bounce as 



Pr 



k^ 
27r2 



487r2M|,, 



(41) 



where Hm = ■\/cr/6 is the absolute value of the Hub- 
ble parameter at the beginning moment of the bouncing 
phase. It should be noted that the Mpi which appears 
in the power spectrum will become rescaled when pmB is 
not equal to 2Afp;cr. 



D. The tensor-to-scalar ratio and the matter 
bounce curvaton scenario 

We can now study the amplitude of tensor perturba- 
tions hij in the f{T) matter bounce. Following (s^], the 
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perturbation equation for the tensor modes can be ex- 
pressed as 

/•• • \ 12HHfTT- 

h,j + mK, - —h,j - —r-i^K, = , (42) 
\ a / t + /,T 

while the tensor modes are transverse and traceless, 
namely 

In the specific case of matter bounce cosmology, equa- 
tion is the same as the standard one since /.tt 
vanishes in the contracting phase. Hence, following the 
analysis done in j46l |. the primordial power spectrum of 
tensor fluctuations is also scale- invariant, but the ampli- 
tude is /i ~ Hm/Mpi, which is of the same order of the 
scalar perturbation. Therefore, the f{T) matter bounce 
scenario suffers from the usual problem of all matter 
bounce models, namely that the ratio of tensor to scalar 
r = Pt/Pc is difficult to accommodate with current ob- 
servations. In particular, from the current CMB data 
[sol ] this ratio is required to be less than 0.2. 

Consequently, in order to make our model consistent 
with cosmological observations on the value of r, it is 
necessary to introduce a mechanism to magnify the am- 
plitude of scalar-type metric perturbations. This issue 
can be resolved by introducing additional light scalar 
fields, as in the matter bounce curvaton scenario [5l| . 
These scalars are able to seed isocurvature fluctuations, 
and then transfer to a scale-invariant spectrum of the 
adiabatic fluctuations during the nonsingular bouncing 
phase, through the so-called kinetic amplification. Thus, 
we obtain a mechanism for enhancing the primordial adi- 
abatic fiuctuations and suppressing the tcnsor-to-scalar 
ratio. 

Following the analysis performed in [5l| we can intro- 
duce one massless scalar field Xj coupling to the back- 
ground matter field through the term g^(j)^x^- Therefore, 
the tensor-to-scalar ratio in our model can be suppressed 
by the kinetic amplification factor in the bouncing phase, 
which can be expressed as r ~ as shown in equation 
(58) of [5l|. Specifically, in our case the kinetic amplifi- 
cation factor T is determined to be ~ 18.66, and so 
r ~ 2.87 X 10~^. Finally, by virtue of the matter bounce 
curvaton mechanism, our model is able to satisfy the con- 
straints from current observations. 



V. CONCLUSIONS AND DISCUSSION 

In this work we investigated the realization of mat- 
ter bounce cosmology in the framework of f(T) grav- 
ity. Considering an explicit scale-factor form, which links 
the matter-like contraction and matter-like expansion 
through a smooth nonsingular bouncing phase, we recon- 
structed the specific form of f{T) that generates it. Our 
analysis has illustrated the possibility of combining the 



/(T) gravity with bounce cosmologies, and thus explor- 
ing an alternative approach to avoid the Big Bang singu- 
larity encountered in the standard infiationary cosmol- 
ogy. Additionally, the constructed /(T) matter bounce 
is free of ghost degrees of freedom and other potential 
problems of many matter bounce scenarios, which is a 
significant advantage. 

At the background level, wc found that the theory of 
torsion gravity is difficult to be verified by experiments. 
This is in agreement with the mechanism described in 
[S^ . where the torsion dynamics is highly degenerated 
with quantum corrections to the classical gravitational 
action caused by vacuum effects. However, going beyond 
the background level, we performed a detailed analysis of 
the evolution of metric perturbations in the present work. 
For our choice of the tetrad wc found that the scalar-type 
metric perturbations posses a single degree of freedom, 
but they have a time-dependent sound speed, a feature 
that is generic in any f{T) scenario. We expect that this 
behavior could be an important signature to be detected 
or constrained by observations. 

The present /(T) matter bounce scenario suffers from 
the usual problem, widely existing in many other bounce 
models, namely it predicts a value for tensor-to-scalar 
ratio larger than observational constraints [i^, [5^ . This 
undesirable feature can be cured by the mechanism of 
matter bounce curvaton, which introduces an additional 
scalar field tracking the background evolution before the 
bounce, but with fluctuations obtaining a kinetic ampli- 
fication during the bouncing phase. 

Finally, we would like to point out that in order to 
establish the existence of matter bounce in f{T) gravity 
and perform a first examination of its basic features, we 
constructed a simplified scenario, without the incorpo- 
ration of the contribution of radiation, dark energy and 
other matter components. Therefore, the investigated 
scenario involves only two parameters, namely the energy 
density at the bouncing point pmB and the bouncing pa- 
rameter a. Through the perturbation analysis we found 
that these two parameters could be tightly constrained 
by CMB observations and experimental bounds on the 
Newton's constant. It would be interesting to construct 
bouncing models in /(T) gravity going beyond this sce- 
nario, in which the contributions of other components 
would be taken into account and various cosmological 
observations would be satisfied. We leave such a con- 
struction of a phenomenologically realistic model and its 
observational constraints for a future investigation. 
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